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We investigate the role of charging effects in a voltage-biased quantum wire. Both the finite range of the Coulomb interaction 
and the long-ranged nature of the Friedel oscillation imply a finite capacitance, leading to a charging energy. While observable 
Coulomb blockade effects are absent for a single impurity, they are crucial if islands are present. For a double barrier, we give 
the resonance condition, fully taking into account the charging of the island. 



PACS numbers: 72.f0.-d, 73.40. Gk 



I. INTRODUCTION 

One-dimensional (ID) quantum wires have attracted 
much interest lately due to theoretical developments and 
advanced fabrication, techniques. By using a special 
split-gate technologytl or cleaved-edge overgrowth] 2 ] high- 
quality ID channels in GaAs-AlGaAs heterostructures 
have recently been fabricated. The typical mean free 
path can be of the order of 10/zm, which brings one close 
to the possibility of experimentally probing the trans- 
port properties of quantum wires, with either none or 
only very few impurities. The most prominent transport 
quantity, the conductance, has been theoretically stud- 
ied in great detaiL, especially for the clean caseBu or for a 
single impurityJj However, little effort has been under- 
taken so far to reveal the nature of charging effectstl in a 
quantum wire. In this paper, we discuss charging effects 
for the simplest case of a spinless single-channel wire. 

The importance of charging effects is established by the 
magnitude of the charging energy e 2 /C. To calculate^his. 
quantity, we utilize standard bosonization methodsETtiJ 
in conjunction with, our recently developed boundary 
condition approacho This formalism allows for a partic- 
ularly simple derivation of the capacitance C of a single 
impurity in the limit of strong impurity backscattering 
strength. It is intuitively clear that the capacitance will 
to a large degree depend on the interactiop-.range R of 
the Coulomb potential in the quantum wireiij Therefore, 
to describe charging effects for realistic experimental se- 
tups characterized by a finite i?, we go beyond the strict 
Luttinger liquid picturetil which has zero range, R = 0, 
and consider an arbitrary screened Coulomb interaction 
potential U(x — y). It turns out that (at least concern- 
ing dc properties) charging effects will only be present if 
the quantum wire contains islands formed by two impu- 
rities. Including the charging contribution, we provide 
the resonance condition for resonant tunneling through a 
double barrier structure. Due to the finite charging en- 
ergy, we feld, a different resonance condition than the one 
predictecB't3 for resonant tunneling in a Luttinger liquid. 



The outline of this paper is as follows. In Sec. |J, 
the boundary condition approach is applied to the case 
of a quantum wire with arbitrary Coulomb interactions 
containing a single impurity. This formalism is used in 
Sec. [II to compute the charge screening cloud around a 
strong impurity, from which one can define the capaci- 
tance C . The resonance condition for resonant tunneling 
through a double barrier is derived in Sec. and the 
additional charging contribution is discussed in detail. 
Finally, some conclusions are offered in Sec. [v} 



II. BOUNDARY CONDITION FORMALISM 

We employ the standard bosonization method, 
which holds in the low-energy regime where only plas- 
mons are well-defined eigenmodes. The electron creation 
operator can then be written in terms of plasmon dis- 
placement fields 9(x) and 4>(x), 



2itvf , 
p=± 



ixp{ipk F x + i*JTr\p8(x) + <f>(x)]} , 



(2.1) 



where oj c — wf^f is the bandwidth (we put % = 1), and 
the sum goes over left- and right-moving pieces (p = ±). 
These fields obey the algebra 

[<P(x),6(y)}=-(i/2) S gn(x-y) , 

such that n = d x (p constitutes the canonical momen- 
tum for the 9 field. The non-interacting ID electron gas 
is then described by the bosonized form of the massless 
Dirac Hamiltonian 



Ho = V -§- |dx[H 2 ( 



x) + (d x 9(x)f 



(2.2) 



Arbitrary screened Coulom b in teractions can be in- 
cluded as follows. From Eq. (2.1), the boson represen- 
tation of the density operator is 



p(x) 



kj,' 

IT 



1 kp _ 

-=d x 6(x) + — cos[2fc F a; + 2y/n9(x)} 
Jit -k 



(2.3) 
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The 2fcF-oscillatory component comes from interference 
of right- asd left-movers and is responsible for a Friedel 
oscillationll-3 in a system with broken translational invari- 
ance. In the standard expression for the electron-electron 
interaction, 



Hj = — J dxdy p(x)U(x 



y)p(y) 



(2.4) 



the interaction among the 2kp parts of and p(y) 
gives rise to electron-electron backscattering JI9 In the fol- 
lowing, we assume that the screened interaction 

U{x - y) = J 7^U k exp[ik(x - y)] 

is sufficiently long-ranged such that Uik F is small. In that 
case, backscattering can be neglected. This is of course 
not an essential assumption but simplifies notation in 
the following. In a homogeneous wire, the interaction 
between the slow and the 2kp part of p averages out due 
to a rapidly oscillating phase factor, i.e., due to momen- 
tum conservation. We then take into account only the 
forward scattering contribution, 



H, 



1 

27 



dxdyd x 9(x) U(x - y) d y 0(y) 



This is t he interaction between the slow components in 
Eq. (p|) . 

Next we discuss how to incorporate coupling of the 
quantum wire to external reservoirs. Incidentally, it has 
become quite customary to use non- interacting ID Lut- 
tinger liquids described by Eq. ( |2.2| ) as a model for the 
external leads BtB There is a serious shortcoming inher- 
ent to such a modelling. If one connects the interact- 
ing wire to ID non-interacting leads, one effectively has 
an inhomogeneous interaction potential U(x,y). In that 
case, the Coulomb interaction between the slow and the 
2fc_p parts of p does not necessarily average to zero since 
momentum is not conserved anymore. In fact, if the in- 
teractions are s witc hed off on a lengthscale 1/kp, we ob- 
tain from Eq. ( |2.4| ) effective potential scatterers at the 
boundaries of the wire of large backscattering strength, 



eff 



1 



- 1 



Hence, from a microscopic point of view, transport would 
practically be suppressed with such leads. This addi- 
tional interaction term has been disregarded in Refs. ^|Jl5[ 
but it is present in a quantum wire connected to ID non- 
interacting leads. 

These difficulties are avoided by the boundary condi- 
tion approach of Ref. |i[ The reservoirs inject currents 
at both ends of the wire. In the spirit of Landauer's ap- 
proach for non-interacting systems the injection of cur- 
rents can be described by Sommerfeld-type radiation con- 
ditions. Coupling to external reservoirs held at a chemi- 
cal potential difference A/x = eU ex leads to the boundary 
conditions 



± 







1 d 



dx vf dt 



(9(x -> Too, t)) = eU cx /2y / Trv F . 

(2.5) 



These conditions have to be imposed at both ends of the 
quantum wire, i.e., far away from all impurities but still 
inside the quantum wire such that bosonization is mean- 
ingful, and for all times t. In general, the effects of an 
applied voltage cannot be captured by adding new terms 
to the Hamiltonian. The boundary condition approach 
holds for arbitrary Coulomb interaction potentials and 
allows for arbitrary arrangements of impurities inside the 
quantum wire. Since we are not directly concerned with 
conductance calculations in the following, we use for sim- 
plicity the T = imaginary-time formalism, where the 
equivalent conditions 



± 







d 



dx vp 8t 



(6(x — > T°o,r = 0)) = eU ex /2-y/nvF 

(2.6) 



are imposed. Here we have picked the time r = by 
convention. 

We first discuss the case of a single impurity located 
at x — 0, and consider the generating functional 

Z(x ,s) = {exp[2^/T:is9(xQ,T = 0)]) . 

One can formally solve for Z by using the auxiliary field 
q(r) = 2y / 7r6>(0, t), where the constraint is enforced by 
a Lagrange multiplier field A(r). The resulting effective 
action is 



S a [6,A,q] 



2 
1 



drdx 



I 



{d T ef + {d x ef 



+ — J dTdxdyd x 6(x)U(x - y)d y 0{y) 
+ V J drcos[g(r)] - 2^is6{x , 0) 
+ i [ dTA(T)[2V^e(0,T)-q(r)} , 



where V denotes the impurity backscattering strength. 
The action is now quadratic in the 9 part, which there- 
fore can be integrated out by solving the classical Euler- 
Lagrange equa tion . This has to be done under the bound- 
ary condition (2.6). 

Analogous to the zero-range Luttinger liquid case 
treated in Ref. ||, this can be achieved by decomposing 
fields into homogeneous and particular parts, 8 = 8h + 9 P 



and A = Ah - 



■Ap. The particular solution < 



Eq. (2.6) and the Euler-Lagrange equation 



— 8%{x, t) + 8%(x,t) 



has to fulfill 



(2.7) 



dy 



U(x - y) 

TTVp 



d 2 J p {y,T) = (2^i/v f )A p {t)5{x) . 
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A p . We make the ansatz 



A solution subject to Eq. ( |2.6| ) requires a r-independent 

e(U ex - if) 



,(x,t) = - 



ep 



2^/ttvf 



it- 



2y^ 

as |x| — > oo for the r-independent 



m - /(o) , 

(2.8) 



where d x f(x) 
function f(x). 

Sw itch ing to Fourier space, one finds from Eqs. (2.7) 
and (Eg) 



eip 
y/nv F 



(1 + Uk/nvp) - k 2 (l + U k /irv F )f k (2.9) 
= (2-y^i /vf) A p . 



The A; = component of Eq. (£Jj) determines the quan- 
tity (p in terms of the zero mode A p of the Lagrange 
multiplier, 



A.,, 



ieip 



Here we have introduced the usual dimensionless Lut- 
tinger liquid parameter g as a-|ineasure of the forward 
scattering interaction strength 

l/g 2 = 1 + U /nv F . 

The non-interacting limit is g = 1, and for repulsive inter- 
actions, one has g < 1. Naturally, a screened interaction 
is characterized by g and the range R (and possibly by 
othe r pa rameters). For R = 0, the k components of 
Eq. (|2.9| ) vanish identically, and /(a;) stays cons tan t. In 
the general finite-range case, we have from Eq. (2.9) the 
simple result 



fk 



apjUo - U k ) 



(2.10) 



which apparently vanishes for the zero-range case where 
Uk = Uq for all k. Here, the plasmon frequency is 



u>k = v F \k\ \/l + Uk/Trv F 



(2.11) 



The homogeneous solution (for U ex = boundary con- 
dition) is easily expressed in terms of the boson propa- 
gators 



F(x, u) — v F I dk 



exp(ikx) 



such that 



9h{x,r) = -- 



Inserting Oh 



2^ 



(2.12) 



e iuT [A h (u>)F(x, u)-sF(x-x , w)] 



back into the action, one observes that 



A/j appears only in quadratic form and can therefore be 



integrated out immediately. After some algebra, one fi- 
nally arrives at 



Z(x,s) = W{x) s2 /eV5K«flp(«.o) 



exp 



2vr v ' F{0,lu) 
(2.13) 

The envelope function W(x) is given by 

duF 2 (x,u})-F 2 (0,u)~ 



W(x) — exp 



2?r 



F(0,w) 



, (2.14) 



and the average over q(ui) = (2tt) 1 J du> q{r) exp(iwr) 
has to be carried out using the action 



du> q(u)q(—cu) eip 



2tt 4F(0,w) 2tt 5 2 
+ V dTCOs[q(r) — ie(U cx — Cd)t 



dr q(r) 



S 



In principle, (p is a fluctuating quantity: One has to 
average over it, because it is the zero- mode of the La- 
grange multipliet-field. Its physical meaning is the four- 
terminal voltage.a In the case of strong impurities consid- 
ered here, Cp is therefore just equal to the two-terminal 
voltage C/ cx . In the following, we restrict ourselves to 
the strong-impurity limit where charging effects are most 
pronounced, and put Cp = U cx . 



III. IMPURITY SCREENING PROFILE AND 
CAPACITANCE 

Let us now analyze the expectation value of the slow 
component of the density, pa(x) — {d x 9)/y/ir. It can 
be obtained by suitable differentiation of the generating 
functional ( |2.13| ), with the result po(x) = d x 9 p (x, 0)/y / 7r, 
since the contributions from the homogeneous solution 
cancel out.EZI Therefore, Eq. (^8|) yields 



Po{x) 



2ttvf 



sgn(x) 



: 9 x f(x) 



(3.1) 



The Fourier transform of f{x) is given in Eq. ( 2.1C ). The 
uniform contribution ~ sgn(x) is due to charging of the 
large shunt capacitances between the quantum wire and 
the metallic (screening and confining) gates, and the cor- 
responding charge on the impurity, Q s — Le 2 U ex / Attvf, 
scales with the total length L of the wire. The observa- 
tion of dc charging effects for a single impurity is re 
dered impossible by this macroscopically large charge 
This can also bp-shown by computing the current-voltage 
characteristics.^ 

The remaining part in Eq. ( |3.l| ) can now be employed 
to provide a microscopic definition of the capacitance C 
of a strong impurity in a quantum wire. Since f(x — > 
±oo) = 0, we have with Q = — e L°° dx[p(x) — p{— x)\/2 
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the simple result Cf r 
itly 

C b = (e/n) 2 



Q/U cx = ef(0)/y/nU ex , or explic- 



dk Up - U k 
2^ w? 



(3.2) 



There is also a contribution C2k F due to the Friedel oscil- 
lation, i.e., the 2Uf component of p in Eq. (|2.3|). While 
the finite-range capacitance Cf r vanishes for the usual 
zero-range Luttinger liquid, Cik F is finite even for R = 
unless one is in the non-interacting limit g = 1. The re- 
sult for C2k F can be found in Ref. ||. In total, since the 
charges simply add up, the capacitance is then given by 
C = + Cik F ■ For the rather long-ranged interactions 
typically present in quantum wires, kpR 3> 1, the slow 
component will dominate, Cf r 3> C2k F ■ 

To give a concrete example, we consider the partic- 
ularly simple form of an exponential interaction with 
dimensionless forward scattering strength u = Uq/ttvf 
such that g = 1/ \/l + u, 



irv F u 

U{x-y) = -^- ex P(" 



y\/R), u k = 



nvpu 



1 



(Rk) 2 ' 
(3.3) 



which allows for explicit analytical calcu latio ns. T he full 
density profile is found from Eqs. (2.1C) and (pO), 



Po(x) 



eU cx 
2t:vf 



sgn(x)[l + ue-W gR } , 



such that the impurity charge at x = is exponentially 
screened on a scale g R. T he finite-range capacitance is 



easily found from Eq. (3.2) as 



e 2 uR 



2'kvf\J\ + u 



(3.4) 



Interestingly, the capacitance is proportional to the in- 
teraction range, Cf r ~ R. This behavior follows imme- 
diately from simple dimensional scaling if the Coulomb 
interaction depends only on \x — y\/R. Furthermore, Cf r 
increases monotonously (but not linearly) with Coulomb 
interaction st rengt h u. 



From Eqs. ( 2.11 ) and ( 2.12 ), one can also evaluat e th e 
boson propagators in closed form for the interaction (3.3). 
The lengthy result can be simplified in two important 
limits. For \u>\ -C Vp/R, the Luttinger liquid result 



gir 

F(x,uj) = -r—r exp[-\gu}x/v F \ 
M 



(3.5) 



is recovered. On the other hand, for |q; | ^ > vf/R, one 
finds the non-interacting result, i.e., Eq. ( ^3.5[ ) with g = 1. 
Apparently there is a new energy scale vf /R associated 
with the interaction range. For frequencies small com- 
pared to this scale, the electrons basically see a zero- 
range (Luttinger liquid) interaction, while for frequencies 
larger than vf/R, the non-interacting behavior is found. 



In that case the electrons are too fast to see each other 
via the finite-range Coulomb interaction. 

Th e F riedel oscillation can be evaluated from 
Eq. ( 2.14 ), since W(x) directly determines the U cx = 
Friedel oscillation for a strong scatterer,E2l 

(p2k F (x)} = -(kF/n)W(x)sm[2k F x\ . 



Splitting up the frequency integration in Eq. (2.14) into 
|w| < vp/R and \u\ > Vp/R, and using the respective 
boson propagators, one finds that approximately 



W{x) ~ (l + 2gx/RY 



2kpx 
l + 2x/R 



For x R, this reproduces the Luttinger liquid re- 
sult, namely an algebraic decay of the Friedel oscilla- 
tion ~ x~ 9 . On the other hand, there is a crossover 
to the faster non-interacting law for x <C R, where the 
Friedel oscillation decays as 1/x. We note that for a 
weak impurity, the situation is more complex since there 
is a competing influence trying to. slow down the Friedel 
oscillation close to the impurity.EJ 

A similar behavior characterizes the conductance. The 
Luttinger liquid power lawsETLl are observed only on en- 
ergy scales small compared to vp/R, while one has the 
non-interacting behavior for larger energy scales. 

Finally, let us comment on the case of an unscreened 
interaction of the formE3 



U(x-y) 



K\J (x - y) 2 + d 2 



2e 2 

K 



(3.6) 



Here, k is the dielectric constant and d denotes the width 
of the wire (kd 1). The capacitance can be ob- 
tainc d fr om Eq. ( |3.2| ). Since the interactions are strong, 
Eq. (3.2) can be simplified to 



C[ r 



71" VF Jl/R 



1 ' d dk 



k 2 



\n[R/d] 
In [fed] 



where a finite interaction range R allows for a controlled 
evaluation of Cf r . The integration yields 



Cfr 



R 



tt 2 v f \n[R/d] 



0{R/ ln 2 [i?/d]) 



(3.7) 



Therefore Cf r diverges ~ Rj lni? as R — * oo for an un- 
screened Coulomb interaction. This underlines the cru- 
cial importance of charging effects in quantum wires with 
a long-ranged l/\x — y\ Coulomb potential.EJ It is note- 
worthy that due to the Rj lni? dependence of the capac- 
itance (3.7), the charge Q s disappears in the unscreened 
case: The interaction range R becomes larger than the 
length of the quantum wire only in the absence of screen- 
ing gates. 
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IV. DOUBLE BARRIER 

The capacitance C and hence Coulomb blockade effects 
can be experimentally observed once islands are present 
in the quant urn r-Wji*e- The simplest case is given by a 
double barrier .BI13I23 We consider two strong impurities 
located at x = id/2, take an infinitesimal two-terminal 
voltage £/ ex and compute the resonance condition as a 
function of the gate voltage (fa coupling to the island 
charge. Naturally, such a set-up might be difficult to re- 
alize experimentally. Nevertheless, the calculations for 
this model show how one can observe charging effects in 
principle. The great merit of such a set-up is that ther- 
modynamic calculations suffice to determine the location 
of the resonances.Eij Since the large barriers confine the 
charge on the island to some integer value n, the reso- 
nance condition at T = is simply E(n) — E(n — 1), 
where E(n) is the energy of the total system with charge 
ne on the island. 

We start from H = Hq + Hi + Hs + Hg, where H s 
describes the impurities and Hq the coupling to the gate 
voltage (pa, 

H G = ep G [6{d/2) - 0(-d/2)]/V5F . 

Since the gate voltage couples to the island charge capac- 
itively, there is no need to resort to a boundary condition 
and one can use the standard term Hq. Furthermore, we 
assume that the impurities are strong enough to pin the 
plasmon displacement fields to the minima of the cosine, 
such that the fields Q{t) and N{f) defined by 



N 



W/2) 
W/2) 



-d/2)]/V5F 
-d/2)]/v^F 



become discrete. Neglecting the Friedel oscillation, the 
charge on the island is n = kpd+N, while Q is associated 
with transport through the island. 

Enforcing the definitions of Q and N by Lagrange mul- 
tiplier fields, one can proceed as before. Inserting the 
solution of the Euler-Lagrange equation into the action 
and integratinffi-put the Lagrange multipliers yields the 
effective actionOj 



it 2 f duj f N{uj)N(-uj) 



S 2 J 2tt \F(Q,oj) 



Q(w)Q(-w) 



-eip G 



F(d,Lu) F(0,u)+F(d,u) 
dTN{r) . 



A = {v F /n 2 ) J dk [1 - cos(fcd)]/^ (4.1) 



= !UL + 9 -ldk[u 



Uk] 



1 — coslkd] 



The first part is the renormalized level spacing derived 
in Refs. ||,13|, while the second part is an additional con- 
tribution due to the charging energy. The condition 
E(N) = E{N — 1) then directly yields the spacing of 
subsequent values of ipQ where one has a resonance, 



e Atp G = 1/A . 



(4.2) 



For interaction range the cos[/cd] in Eq. (4.1) 

does not contribute, and one re cov ers the single-impurity 
capacitance Cf r defined in Eq. ( |j.2| ) . Taking into account 
also the Friedel oscillation contribution Cik F , on e has the 
resonance condition in a particularly simple form, 



eg 2 A(p G = 



d 

TTVf 



2C 

~Z2 



(4.3) 



Remarkably, the capacitance will decrease the spacing 
Aipc and therefore AffH is smaller than predicted for 
the zero-range model.ElEj At first sight, this might ap- 
pear counter-intuitive because charging effects suppos- 
edly increase the spacing.LJ However, since A i s al ways 
diminuished by repulsive interactions, see Eq. (4.1), in- 
teractions per se will alw ays enhance the spacing Atpo- 
The factor g 2 in Eq. (4.3) tends to increase A<pc while a 
finite capacitance C decreases the spacing again. In total, 
however, one is still left with an enhanced spacing com- 
pared to the resonant-tunneling value eAipG — irvp/d 
of a non-interacting wire. As a simple example, we 
consider Eq . (4.3) for the exponential i nter action po- 
tential ([O) leading to the capacitance ( |3.4| ). In that 
case, Eq. ( |4.3| ) becomes with the dimensionless interac- 
tion strength u = Uq/ttvf 



eAipG 



TTVp(l + u) 

d + Ru/^/TT 



For a given interaction range R, increasing the interac- 
tion strength u will always lead to a larger spacing eAtpc 
compared to the value irvp/d arising from the bare level 
spacing of the island. On the other hand, for given u, 
increasing R decreases Aipo f° r the reasons discussed 
above. 



The mode TV is gapped while Q is ungapped and hence 
irrelevant with respect to charging effects. The low- 
frequency sector can therefore be described by an effec- 
tive energy 

E{N) = N 2 /2A + etp G N , 

where we have introduced the quantity A — [F(0, 0) — 
F(d, 0)]/7r 2 whose explicit form is 



V. CONCLUSIONS 

In this paper, we have employed the boundary con- 
dition formalism to investigate charging effects in one- 
dimensional quantum wires. Here, charging effects have 
two distinct origins: (i) The finite range of the screened 
Coulomb interaction potential, and (ii) the long-ranged 
nature of the Friedel oscillation in a Luttinger liquid. The 
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full density profile around an impurity in the presence 
of external voltage sources can be computed using the 
bosonization method under appropiate boundary condi- 
tions. Using that result, one can infer the value of the 
capacitance of the impurity. While this capacitance and 
hence charging effects do not seem to have consequences 
for dc transport properties in the case of a single impu- 
rity, the condition for resonant tunneling through a dou- 
ble barrier is modified compared to previous estimates if 
charging is properly taken into account. 
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